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Temperature  has  a  huge  influence  on  supercapacitor  cells  and  modules  ageing.  Consequently,  thermal 
management  is  a  key  issue  concerning  lifetime  and  performance  of  supercapacitor  modules.  This  paper 
presents  thermal  modeling  and  heat  management  of  supercapacitor  modules  for  vehicle  applications.  The 
thermal  model  developed  is  based  on  thermal-electric  analogy  and  allows  the  determination  of  super¬ 
capacitor  temperature.  Relying  on  this  model,  heat  management  in  supercapacitor  modules  was  studied 
for  vehicle  applications.  Thus,  the  modules  were  submitted  to  real  life  driving  cycles  and  the  evolution  of 
temperatures  of  supercapacitors  was  estimated  according  to  electrical  demands.  The  simulation  results 
show  that  the  hotspot  is  located  in  the  middle  of  supercapacitors  module  and  that  a  forced  airflow  cooling 
system  is  necessary. 
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1.  Introduction 

Heat  production  in  supercapacitor  is  related  exclusively  to  Joule 
losses.  The  supercapacitors  support  currents  up  to  400  A  or  more 
depending  on  cell  capacitance  and  used  technology.  The  repet¬ 
itive  charge  and  discharge  cycles  of  the  supercapacitor  cause  a 
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significant  warming  even  though  the  equivalent  series  resistance 
value  is  around  the  m£2  according  to  the  capacitance.  Several 
authors  showed  that  the  supercapacitor  ESR  varies  according 
to  the  temperature  [1-3].  In  [4]  the  authors  have  studied  the 
effect  of  the  temperature  and  the  voltage  on  the  supercapacitors 
ageing.  They  have  established  a  model  which  allows  analyzing  self- 
accelerating  degradation  effects  caused  by  elevated  voltages  and 
temperatures.  In  the  reference  [5]  the  authors  have  studied  and 
modeled  the  temperature  effect  on  the  supercapacitor  self  dis¬ 
charge. 


0378-7753 /$  -  see  front  matter  ©  2009  Elsevier  B.V.  All  rights  reserved. 
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Nomenclature 

cp  specific  heat  capacity  (J  kg-1  °C_1 ) 

C  supercapacitor  capacitance  (F) 

Cth  thermal  capacitance  (J°  C-1 ) 

Cth_ i  thermal  capacitance  of  the  layer  number  i  (J°  C-1 ) 

Qnsuiating ^ -layer  thermal  capacitance  of  the  insulating  layer 

(re-1) 

Qnetai  -case  thermal  capacitance  of  the  metal  case  (J°  C-1 ) 

dr  elementary  thickness  (m) 

dRth  elementary  thermal  resistance  (£2) 

EIS  electrochemical  impedance  spectroscopy 

ESR  supercapacitor  equivalent  series  resistance  (£2) 

ES  energy  storage 

EVs  electric  vehicles 

/  frequency  (Hz) 

FCVs  fuel  cell  vehicles 

G  electrical  conductance  ( ) 

Gth  thermal  conductance  Tme a 

h  convection  heat  transfer  coefficient  (W m-2  °C_1 ) 

HyHEELS  hybrid  high  energy  electrical  storage 
HEVs  hybrid  electric  vehicles 

I  current  flux  (A) 

/( t)  RMS  current  value  (A) 

Im(z)  imaginary  component  of  the  supercapacitor 

impedance  response 
J  current  flux  density  (A  irr 2 ) 

L  hollow  cylinder  length  (m) 

n  number  of  layers  of  supercapacitor 

P(t)  total  power  dissipated  by  Joule  effect  (W) 

R  hollow  cylinder  inner  radius  (m) 

^electric  electrical  resistance  (£2) 

Rs  hollow  cylinder  outer  radius  (m) 

Re(z)  real  component  of  the  supercapacitor  impedance 

response 

^convection  convection  heat  transfer  resistance  (£2)  or 

^cw-1) 

Rth  thermal  resistance  (°C  W_1 ) 

Rth_i  thermal  resistance  of  the  layer  number  i  (°C  W-1 ) 

insulating -layer  thermal  resistance  of  the  insulating  layer 

(^CW-1) 

^metal-case  thermal  resistance  of  the  metal  case  (°C  W_1 ) 

Ssc  heat  exchange  surface  of  the  supercapacitor  (m2) 

T  hollow  cylinder  inner  surface  temperature  (°C) 

Ta  ambient  temperature  (°C) 

^Hi-mea  measured  temperatures  (°C) 

THi-sin  simulated  temperatures  (°C) 

Tinea  measured  temperature  (°C) 

Tmin  minimum  temperature  (°C) 

Tmax  maximum  temperature  (°C) 

Ts  hollow  cylinder  outer  surface  temperature  (°C) 

Tsin  simulated  temperature  (°C) 

AT  temperature  difference  (°C) 

AV  potential  difference  (V) 

Greek  letters 
0  heat  flux  (W) 

cp  heat  flux  density  ( W  irr 2 ) 

A  thermal  conductivity  (W m-1  °C_1 ) 

fi  density  (Kg  m-3) 

a  electrical  conductivity  THi-Sin 

r  system  time  constant  (s) 


0.01  0.1  1  10  100  1000 

Frequency  (Hz) 


Fig.  1.  BCAP1500F  and  BCAP310F  capacitance  as  function  of  frequency  with  a  bias 
voltage  respectively  of  2.7  V  and  2.5  V  and  a  temperature  of  20  °C. 

This  rise  in  temperature  can  have  the  following  consequences: 

•  The  deterioration  of  the  supercapacitor  characteristics,  especially 
ESR,  self  discharge  and  lifetime  [4,5],  which  affect  its  reliability 
and  its  electrical  performance. 

•  The  pressure  inside  the  supercapacitor  is  increased. 

•  A  premature  aging  of  metal  contacts,  in  fact  the  repetitive  heating 
and  significant  temperatures  can  deteriorate  rapidly  the  terminal 
connections  of  the  supercapacitor. 

•  The  evaporation  of  the  electrolyte  and  hence  the  destruction  of 
the  supercapacitor  if  the  temperature  exceeds  81.6  °C  the  boiling 
point  of  the  electrolyte. 

Therefore,  it  is  important  to  know  and  understand  the  heat 
behavior  of  supercapacitor  cells  and  modules.  This  leads  to  an  esti¬ 
mation  of  the  space-time  evolution  of  the  temperature. 

This  paper  deals  with  the  thermal  modeling  of  supercapacitors 
and  heat  management  in  supercapacitor  modules.  The  originality 
of  this  work  is  based  on  the  integration  of  thermocouples  located 
inside  the  supercapacitor  during  its  manufacturing  by  Maxwell 
technologies.  Cooling  systems  were  also  studied  for  supercapacitor 
modules  subjected  to  a  driving  cycle. 

2.  Supercapacitor  electric  characterization 

The  Maxwell  BCAP310F  and  BCAP1500F  supercapacitors  used 
in  this  study  are  based  on  activated  carbon  technology  and 
organic  electrolyte.  These  devices  were  characterized  using 
the  electrochemical  impedance  spectroscopy  (EIS).  This  allows 
the  determination  of  the  supercapacitor  real  and  imaginary 
components  of  the  impedance  response.  It  assumed  that  the  super¬ 
capacitor  capacitance  C  and  the  series  resistance  (ESR)  are  deduced 
from  the  experimental  results  respectively: 


2jTlm{z)f 
ESR  =  Re{z) 

where  Im(z)  and  Re(z)  are  respectively  the  imaginary  part  and  the 
real  part  of  the  supercapacitor  impedance  and /is  the  frequency. 

Figs.  1  and  2  represent  the  BCAP310F  and  the  BCAP1500F  capac¬ 
itance  and  ESR  as  a  function  of  frequency.  At  low  frequency, 
the  capacitance  is  maximum,  for  example  at  10  mHz  the  capac¬ 
itance  value  is  in  order  of  1660  F  for  the  BCAP1500F  and  315  F 
for  the  BCAP310F.  At  50  mHz  the  ESR  value  is  in  order  of  1 
for  BCAP1500F  and  5.2  m^  for  BCAP310F.  The  BCAP310F  ESR  is 
relatively  high  because  this  device  was  fabricated,  by  Maxwell  tech¬ 
nologies,  especially  for  these  thermal  tests;  it  is  including  four 
thermocouples  type  K  inside. 


(1) 

(2) 
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Fig.  2.  BCAP1500F  and  BCAP310F  series  resistance  as  function  of  frequency  with  a 
bias  voltage  respectively  of  2.7  V  and  2.5  V  and  a  temperature  of  20  °C. 
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Fig.  3.  BCAP310F  equivalent  series  resistance  as  function  of  temperature. 

The  BCAP1500F  capacitance  and  ESR  were  calculated  in  DC  state. 
ESR  and  C  are  respectively  about  1.07  mC2  and  1525  F. 

For  supercapacitor  thermal  behavior,  the  device  was  character¬ 
ized  by  using  the  EIS  for  different  temperature.  Fig.  3  represents  the 
Maxwell  BCAP0310F  ESR  variations  according  to  the  temperature. 
The  ESR  increases  at  negative  temperature  values.  The  ESR  variation 
is  higher  for  negative  temperature  than  for  positive  one.  This  is  due 
to  the  fact  that  the  electrolyte’s  ionic  resistance  is  strongly  influ¬ 
enced  by  the  temperature.  Above  0°C  ESR  varies  slowly  with  the 
temperature.  Below  0°C  the  temperature  dependence  is  stronger. 
Higher  ESR  is  due  to  the  increase  of  the  electrolyte’s  viscosity  at 
low  temperatures  limiting  ionic  transport  speed  which  increases 
the  resistance  of  the  electrolyte. 

In  the  case  of  the  capacitance,  the  experimental  results  show 
that  the  capacitance  is  lower  at  negative  temperature  as  shown 
in  Fig.  4.  For  example,  at /=  10  mHz  there  is  no  variation  of  the 


Temperature  (°C) 


Fig.  4.  BCAP310F  capacitance  evolution  according  to  the  temperature  for  10  mHz 
and  lOOmHz. 


capacitance  with  temperature.  At  lOOmHz,  C=335F  at  -20 °C 
whereas  C=361  F  at  20  °C.  At  negative  temperature,  the  superca¬ 
pacitor  capacitance  decreases  with  temperature. 

In  conclusion,  it  is  clear  that  the  supercapacitor  electric  perfor¬ 
mances  and  lifetime  depend  on  the  temperature. 

3.  Thermal  modeling  of  supercapacitors 

In  order  to  establish  a  thermal  model  of  a  supercapacitor  cell,  it 
is  important  to  expose  its  geometric  structure.  The  structure  of  a 
basic  cell  is  cylindrical.  The  technology  achievement  is  identical  to 
that  used  for  conventional  capacitors.  Thus,  there  is  a  basic  multi¬ 
layer  that  consists  of  a  positive  electrode  and  a  negative  electrode 
of  activated  carbon.  They  are  electrically  isolated  with  a  separa¬ 
tor  placed  between  them.  The  ensemble  forms  a  layer  which  is 
rolled  several  times  and  then  placed  in  a  metal  case  and  impreg¬ 
nated  with  an  organic  electrolyte.  The  two  electrodes  are  metalized 
and  connected  to  the  outside  (+)  and  (-)  terminal  connections  of 
supercapacitor. 

Refs.  [6-9]  present  thermal  modeling  of  supercapacitors.  Differ¬ 
ent  methods  are  used,  heat  equation  resolution,  electric-thermal 
analogy.  The  models  developed  allow  the  determination  of  super¬ 
capacitor  temperature. 

As  it  described  before,  the  structure  is  cylindrical  geometry.  The¬ 
oretical  study  and  experimental  results  show  that  the  temperatures 
of  a  layer  are  uniform;  however  the  temperature  of  each  layer  is  dif¬ 
ferent.  Thus,  the  thermal  modeling  can  be  studied  by  calculating  the 
total  thermal  resistance  and  capacitance  of  the  supercapacitor. 

3.1.  Theoretical  recall 

Consider  a  hollow  cylinder  with  an  inner  radius  R,  an  outer 
radius  Rs,  a  length  L,  thermal  conductivity  A,  specific  heat  capacity 
cp  and  density  pi.  The  inner  surface  is  maintained  at  a  temperature 
T  and  the  outer  surface  at  temperature  Ts. 

The  elementary  thermal  resistance  and  thermal  capacitance  of 
an  elementary  thickness  dr  is  expressed  by: 


dRtii  =  ]  d' 

A  2irrL 

(3) 

dCth  =  pLcp2nLrdr 

(4) 

Therefore,  the  total  resistance  and  capacitance  are  given  by: 

Rs 

1  f dr  1  fRs\ 

Kth  -  2 ttAL  J  r  ~  2 jtXL  n  V  R  ) 

R 

(5) 

Cth  =  licp7zL{Rs2  -  R2) 

(6) 

3.2.  Case  of  a  supercapacitor 

In  order  to  calculate  the  total  thermal  resistance  Rth  and  thermal 
capacitance  Cth  of  the  supercapacitor,  each  layer  in  the  superca¬ 
pacitor  is  modeled  by  its  thermal  resistance  Rth_i  and  its  thermal 
capacitance  Cth_i.  Rth  and  Cth  are  given  by: 

n 

^th  —  ^insulatingJayer  +  ^metaLcase  T  ^  ^th_i 

i=l 

(7) 

n 

Qh  —  QnsulatingJayer  +  OnetaLcase  ^  ^Qh_i 

i=l 

(8) 

Where:  Rmetal-case*  ^insulating -layer*  Qnetal-case  <*nd  Qnsulating -layer 

are  the  thermal  resistance  and  thermal  capacitance  of  the  metal 
case  and  insulating  layer  and  n  is  the  number  of  layers. 
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Table  1 

Thermal-electrical  equivalences. 

Thermal  parameters 

Electrical  parameters 

Thermal  conductivity  (A)  (W m_1  °C_1 ) 
Temperature  difference  (AT)  (°C) 

Heat  flux  density  ( (p )  (W  m-2 ) 

Thermal  resistance  ( Rth )  (0CW_1) 
Thermal  conductance  (Gth)  (W°C-1) 
Heat  flux  (0)  (W) 

Thermal  capacitance  (Cth)  (J°C-1 ) 

Electrical  conductivity  ( a  )  (£2_1  m-1 ) 
Potential  difference  (AV)  (V) 

Current  flux  density  (])  (A  m-2 ) 
Electrical  resistance  (^electric)  (^) 
Electrical  conductance  (G)  (£2_1) 
Current  flux  (/)  (A) 

Electrical  capacitance  (C)  (F) 

Based  on  this  analysis,  a  Matlab/Simulink®  simulation  model 
was  developed  in  order  to  calculate  the  Rt h  and  Cth  of  a  supercapac¬ 
itor  cell.  Calculated  values  were  compared  to  experimental  values 
and  the  simulation  model  was  validated.  This  application  permits 
to  calculate  the  evolution  of  the  temperature  in  each  layer  of  the 
supercapacitor  cell. 

The  thermal  modeling  of  a  supercapacitor  is  done  using  the 
thermal-electric  analogy.  Table  1  shows  the  equivalences  between 
thermal  and  electrical  parameters. 

Fig.  5  shows  the  thermal-electric  model  of  the  supercapacitor. 
The  conduction  heat  transfer  mode  is  the  most  important  within  a 
supercapacitor.  The  other  two  modes  of  heat  transfer  (radiation  and 
convection)  are  neglected.  However  the  convection  heat  transfer 
mode  between  ambient  air  and  the  outer  surface  of  the  supercapac¬ 
itor  is  taken  into  account.  The  thermal  model  gives  the  evaluation  of 
the  temperature  of  the  supercapacitor  depending  on  the  electrical 
power,  the  ambient  temperature  and  the  convective  heat  transfer 
coefficient. 

The  total  power  dissipated  in  the  supercapacitor  is  given  by: 
P(t)  =  ESR  x  /2(t)  (9) 

Where:  ESR  is  the  equivalent  series  resistance  of  the  superca¬ 
pacitor  and  J(t)  is  the  RMS  current  value  passing  through  it. 

The  resistance  RConvection  represents  the  heat  transfer  between 
the  surface  of  the  supercapacitor  and  the  ambient  air.  Its  value 
depends  on  the  convective  heat  transfer  coefficient  h  and  the  heat 
exchange  surface  of  the  supercapacitor  Ssc.  This  coefficient  can  be 
calculated  by  determining  the  system  time  constant  r. 

Where  .  X  =  ^convection  x  Qh  ^convection  —  (10) 

^convection  can  also  be  calculated  by: 

^convection  —  (H) 


Table  2 

Parameters  values  of  310  F  and  1500  F  supercapacitor  cells. 


Parameters 

310  F 

1500  F 

Ta 

22.5  (°C) 

17.5  (°C) 

^convenction 

11  (°C/W) 

3.5  (°C/W) 

^th 

6.5  (°C/W) 

4.2  (°C/W) 

Cth 

44  (J/°C) 

268(J/°C) 

supercapacitor  cell,  with  four  thermocouples  placed  inside  was 
developed  in  order  to  validate  the  thermal  model.  The  model  was 
also  validated  for  a  1500  F  cell. 

Table  2  shows  the  parameters  values  used  in  the  model  for  a 
310  F  and  1500  F  supercapacitor  cells. 

As  explained  previously,  ^convection  depends  on  the  convec¬ 
tive  heat  transfer  coefficient  h  (here  natural  convection)  and  the 
heat  exchange  surface  of  the  supercapacitor  5SC.  It  is  normal  that 
^convection  is  smaller  in  case  of  the  1500  F  cell  because  RConvection 
and  the  heat  exchange  surface  are  inversely  proportional,  and  the 
outer  surface  of  the  1500  F  is  bigger  than  the  310  F  one. 

In  order  to  study  the  condition  of  repetitive  and  extreme  stress, 
repetitive  charge  and  discharge  cycles  between  (1.25  V  and  2.5  V) 
and  (1.35  V  and  2.7  V)  respectively  with  30  A  and  75  A  constant  cur¬ 
rents  are  applied  to  the  310  F  and  1500  F  cells,  respectively.  The 
supercapacitors  were  suspended  in  the  air  to  avoid  any  contact 
with  other  materials.  The  measurement  is  carried  out  with  free 
convection. 

Fig.  6  shows  the  measured  temperatures  THi_mea  (i  =  h  2,  3,  4) 
and  those  obtained  by  simulation  THi_sim  f°r  the  310  F  cell.  Tm ,  TH 2, 
TH3  and  TH 4  represent  the  temperatures  of  the  four  thermocouples 
placed  inside  the  BCAP310F  supercapacitor.  The  four  Thermocou¬ 
ples  are  placed  at  the  middle  of  the  axial  position.  In  the  radial 
position,  the  first  thermocouple  which  is  represented  by  the  tem¬ 
perature  THi  is  placed  after  three  turns  of  jelly  roll,  the  second 
thermocouple  IH 2  is  placed  after  six  turns  of  jelly  roll,  the  third 
thermocouple  TH3  is  placed  after  nine  turns  of  jelly  roll  and  the 
fourth  thermocouple  TH 4  is  placed  after  12  turns  of  jelly  roll.  The 
first  phase  represents  the  warming  phase  of  the  cell  which  is  sub¬ 
jected  to  the  charge  and  discharge  cycles  during  around  40  min. 
Then,  during  the  second  phase  the  process  is  stopped  and  the  tem¬ 
perature  of  the  supercapacitor  is  decreased  by  natural  convection 
to  reach  the  ambient  temperature  which  is  about  22.5  °C. 

Fig.  7  shows  the  evolution  of  the  outside  surface  temperature 
of  the  1500  F  supercapacitor.  The  warming  phase  is  about  133  min 


For  an  optimal  design  of  the  supercapacitor  and  its  cooling  sys¬ 
tem  in  case  of  supercapacitor  modules,  it  is  advantageous  to  know 
the  spatial  and  temporal  evolution  of  temperature.  A  special  310  F 


Ta  R  convection 


P  B 


AM — , 
AM— 


-t> 


Fig.  5.  Thermal-electric  model  of  the  supercapacitor. 


Fig.  6.  Evolution  of  measured  and  simulated  temperatures  of  310  F  cell  as  function 
of  time.  Tm,  Tm,  Tm  and  Tm  represent  the  temperatures  of  the  four  thermocouples 
placed  inside  the  BCAP310F  supercapacitor.  The  four  Thermocouples  are  placed  at 
the  middle  of  the  axial  position.  In  the  radial  position,  Tm  is  placed  after  3  turns  of 
jelly  roll,  Tm  is  placed  after  6  turns  of  jelly  roll,  TH  3  is  placed  after  9  turns  of  jelly  roll 
and  TH4  is  placed  after  12  turns  of  jelly  roll. . 
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Fig.  7.  Evolution  of  measured  and  simulated  temperatures  of  1500  F  cell  as  function 
of  time. 

where  the  supercapacitor  is  charged  and  discharged  at  75  A  as  con¬ 
stant  current.  The  ambient  temperature  is  around  17.5  °C. 

These  results  show  a  good  correlation  between  experimentation 
and  simulation.  Good  agreements  were  obtained  with  10  A  and  20  A 
as  a  current  value  of  charge  and  discharge  cycles  as  well. 

4.  Thermal  management  of  supercapacitor  modules 

Supercapacitor  thermal  management  is  important  in  order  to 
optimize  the  performance  and  to  reduce  the  life-cycle  costs  of 
advanced  fuel  cell,  hybrid  electric  and  electric  vehicles  (FCVs,  HEVs, 
and  EVs).  Both,  temperature  and  temperature  uniformity  signifi¬ 
cantly  affect  the  performance  and  life  of  energy  storage  (ES)  devices 
and  vehicles.  That  is  why  the  European  project  ElyHEELS  (hybrid 
high  energy  electrical  storage)  conducts  research  and  development 
in  the  thermal  management  supercapacitors  to  optimize  perfor¬ 
mance  and  extend  life  of  the  supercapacitors,  the  fuel  cell  and  the 
whole  power  train  of  the  vehicle.  In  fact,  ElyEIEELS  provides  a  super¬ 
capacitor  energy  storage  system  for  the  use  in  hybrid  and  fuel  cell 
vehicles,  which  satisfies  all  properties  necessary  to  make  it  an  inte¬ 
grative  component. 

In  order  to  study  the  thermal  management  of  supercapacitor 
modules,  it  is  interesting  to  study  the  influence  of  the  configura¬ 
tion  of  the  supercapacitor  cells  in  the  module.  A  simulation  model 
in  Matlab/Simulink®  is  developed  to  study  by  simulation  the  ther¬ 
mal  management  in  supercapacitor  modules.  Each  cell  is  modeled 
by  a  thermal-electric  model.  The  convection  heat  transfer  mode 
between  cells  is  modeled  by  resistances.  The  convection  from  a  cell 


Table  3 

Characteristics  of  Maxwell  1500  F  supercapacitor. 


Parameters 

Provided  by  Maxwell 

Measured 

Operating  temperature  range 

-40  °C,  +65  °C 

- 

Equivalent  series  resistance 

ESR  =  0.47  (m£2) 

ESR  =  0.39(mft) 

Thermal  resistance 

Rth  =  4.5(°C/W) 

Rth  =  4.2(°C/W) 

Thermal  capacitance 

- 

Cth  =  268  (J/°C) 

is  modeled  by  resistances  which  depend  on  the  heat  exchange  sur¬ 
face  and  the  convective  heat  transfer  coefficient  (the  number  of 
resistances  is  depends  on  the  number  of  the  neighboring  cells).  The 
convection  between  cells  and  air  is  also  modeled  by  a  voltage  source 
which  represents  the  ambient  temperature,  and  a  resistance  which 
represents  the  convection  heat  transfer.  The  modeling  takes  into 
account  the  convective  heat  transfer  coefficient  which  changes  if 
the  convection  is  natural  or  forced.  Fig.  8  gives  an  idea  about  the 
thermal  modeling  of  a  six  supercapacitors  module  and  the  heat 
transfer  flows  between  the  air  and  the  cells  and  between  the  cells 
themselves. 

This  simulation  model  makes  it  possible  to  dimension  the  super¬ 
capacitor  cooling  system,  if  it  is  necessary.  This  is  in  order  to 
maintain  the  temperature  of  the  module  in  the  operating  temper¬ 
ature  range  given  by  the  manufacturer. 

4.1.  Simulation  specifications 

The  supercapacitor  used  in  the  simulation  is  a  BCAP  1500  Farad 
(Maxwell  technology);  its  characteristics  are  shown  in  Table  3: 

The  module  studied  in  this  paper  is  composed  of  20  cells  of 
1500  F.  The  cells  are  arranged  into  four  rows  of  five  supercapac¬ 
itors  each.  The  simulation  was  carried  out  using  the  simulation 
model  presented  previously.  Each  supercapacitor  is  represented  by 
its  thermal  model.  The  manufacturer  recommends  an  operating 
temperature  range  between  -40  and  65  °C.  The  study  of  thermal 
management  gives  the  evolution  of  the  temperatures  in  different 
points  of  the  module  and  determines  if  the  module  needs  a  cooling 
system.  The  simulation  model  permits  also  to  size  the  cooling  sys¬ 
tem  when  it  is  necessary.  The  simulation  takes  into  account  natural 
or  forced  convection  and  conduction  heat  transfer  modes. 

The  application  permits  to  subject  the  module  to  any  power  pro¬ 
file.  For  vehicle  and  transport  application,  it  is  interesting  to  study 
the  evolution  of  the  temperatures  of  supercapacitor  modules  when 
the  module  is  subjected  to  a  real  life  driving  cycle.  In  this  study  the 
Scania  and  BMW  driving  cycles,  respectively  for  heavy  duty  vehicle 
and  for  passenger  car  applications,  are  considered.  The  character¬ 
istics  of  each  cycle  are  presented  below. 
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Fig.  8.  Thermal-electric  model  of  six  supercapacitors  module. 
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Fig.  9.  Supercapacitors  modules’  current  for  Scania  driving  cycle  as  a  function  of 
time. 


4.1.1.  Scania  heavy  duty  vehicle  driving  cycle  characteristics: 

The  total  time  of  the  cycle  is  22.766  min  (1366  s).  Fig.  9  shows 
the  current  of  the  supercapacitor  modules  during  a  cycle.  The  power 
losses  in  the  supercapacitor  modules  depend  on  this  current  profile. 

4.1.2.  BMW  passenger  car  driving  cycle  characteristics: 

The  total  time  of  the  cycle  is  11.975  min  (718.5  s).  Fig.  10  shows 
the  current  of  the  supercapacitor  modules  during  a  cycle.  The  Joule 
effect  in  the  supercapacitor  modules  is  caused  by  this  current  pro¬ 
file. 

4.2.  Simulation  results 

The  simulations  are  carried  out  for  two  different  initial  ambi¬ 
ent  temperatures  which  are  25  °C  and  50  °C  respectively,  and  for 
one  and  more  cycles  in  order  to  get  the  steady  states  (Scania: 
30  cycles  =  11.383  h,  BMW:  60  cycles  =  11.975  h)  with  natural  and 
assisted  convection. 

4.2.1.  Scania  heavy  duty  vehicle  driving  cycle 

For  one  driving  cycle,  supercapacitors  module  does  not  require 
a  cooling  system.  However,  in  practice  supercapacitors  module  is 
used  usually  for  more  than  one  driving  cycle. 

Figs.  11  and  12  show  the  evolution  of  the  maximum  temper¬ 
ature  of  the  supercapacitors  module  as  a  function  of  time  with 
natural  and  assisted  convection  and  with  25  °C  and  50  °C  as  ambi¬ 
ent  temperatures,  respectively.  The  maximum  temperature  of  the 
supercapacitors  module  is  located  in  the  middle  of  the  module. 
In  natural  convection,  the  temperature  is  higher  than  65  °C,  which 
means  that  the  supercapacitors  module  requires  a  cooling  system. 


Time  [s] 


Fig.  10.  Supercapacitors  modules’  current  for  the  BMW  driving  cycle  as  a  function 
of  time. 


Fig.  11.  Evolution  of  the  maximum  temperature  in  the  module  as  function  of  time 
(30  cycles,  ambient  temperature  Ta  =  25  °C,  natural  and  forced  convection)  with  Sca¬ 
nia  driving  cycle. 


Fig.  12.  Evolution  of  the  maximum  temperature  in  the  module  as  function  of  time 
(30  cycles,  ambient  temperature  Ta  =  50  °C,  natural  and  forced  convection)  with  Sca¬ 
nia  driving  cycle. 


This  last  was  introduced  in  the  simulation  by  changing  the  convec¬ 
tive  heat  transfer  coefficient.  In  forced  convection,  it  is  clear  that 
the  temperatures  are  less  than  65  °C. 

4.2.2.  BMW  passenger  car  driving  cycle 

Fig.  13  shows  the  evolution  of  the  maximum  and  minimum  tem¬ 
peratures  of  the  supercapacitors  module  as  a  function  of  time  for 


Time  [mn] 

Fig.  13.  Evolution  of  the  minimum  and  maximum  temperatures  in  the  module  as 
function  of  time  (60  cycles,  ambient  temperature  Ta  =  25  °C,  natural  convection)  with 
BMW  driving  cycle. 
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Fig.  14.  Evolution  of  the  maximum  temperature  in  the  module  as  function  of  time 
(60  cycles,  ambient  temperature  Ta  =  50°C,  natural  and  forced  convection)  with 
BMW  driving  cycle. 


Table  4 

Maximum  and  minimum  temperatures  in  all  cases. 


Case 

1  Cycle 

30  Cycles 

Tmin  (°C) 

7m ax  (°C) 

Tmin  (°C) 

Tmax  (°C) 

Tmax  >65  (°C) 

Scania  cycle 

Ta  =  25  (°C) 

29.21 

32.21 

39.42 

72.41 

After  10  cycles 

Ta  =  50  (°C) 

54.21 

57.21 

64.42 

97.4 

After  2  cycles 

Ta  =  25  (°C)  + Cooling 

29.42 

39.52 

Never 

Ta  =  50  (°C)  + Cooling 

54.42 

64.52 

Never 

1  Cycle 

60  Cycles 

Tmin  (°C) 

Tmax  (°C) 

Tmin  (°C) 

Tmax  (°C) 

Tmin  >  65  (°C) 

BMW  cycle 

Ta  =  25  (°C) 

26.7 

27.75 

35 

58.5 

Never 

Ta  =  50  (°C) 

51.7 

52.75 

60 

83.5 

After  7  cycles 

Ta  =  50  (°C)  +  Cooling 

53 

60.25 

Never 

The  bold  values  mean  that  the  temperature  is  limit  or  higher  than  65  °C. 


60  cycles  with  natural  convection  and  25  °C  as  an  ambient  temper¬ 
ature.  Fig.  14  shows  the  evolution  of  the  maximum  temperature  of 
the  supercapacitors  module  as  a  function  of  time  with  natural  and 
assisted  convection  and  with  50  °C  as  ambient  temperature.  In  nat¬ 
ural  convection,  the  maximum  temperature  is  higher  than  65  °C;  a 
cooling  system  was  required,  sized  and  verified  as  shown  in  Fig.  14. 

4.2.3.  Required  cooling  system: 

In  this  application,  the  convective  heat  transfer  coefficient  h  in 
natural  convection  is  approximately  12  W  rrr2  °C_1 .  For  assisted  or 
forced  convection  a  higher  convective  coefficient  is  required  which 
is  about  three  times  more  than  the  natural  one.  This  is  equivalent 
to  an  air  speed  of  2.5  m  s-1  approximately. 

Table  4  recapitulates  all  the  cases.  In  conclusion,  for  BWM  driv¬ 
ing  cycle  the  supercapacitors  module  does  not  require  a  cooling 


system,  if  the  ambient  temperature  is  25  °C.  Flowever,  the  cool¬ 
ing  system  is  necessary  when  the  ambient  temperature  is  in  order 
of  50  °C.  For  Scania  heavy  duty  vehicle  cycle,  a  cooling  system  is 
required  because  the  module  temperature  exceeds  the  maximum 
temperature  given  by  the  manufacturer  (65  °C). 

5.  Conclusion 

In  this  paper,  a  supercapacitor  thermal  model  is  developed. 
This  last  is  based  on  the  thermal-electric  analogy.  Thermal 
resistance  and  thermal  capacitance  have  been  calculated  for  a 
cylindrical  geometry.  The  model  is  validated  using  experimental 
results  of  the  BCAP310F.  This  last  is  a  special  supercapaci¬ 
tor  cell  including  four  thermocouples  inside  manufactured  by 
Maxwell  technologies.  The  thermal  model  was  adapted  for 
BCAP1500F.  Simulation  and  experimental  results  for  BCAP1500F 
and  BCAP310F  are  in  good  agreement.  The  thermal  model 
developed  in  this  paper  was  used  for  heat  management  of  super¬ 
capacitor  modules  for  Scania  heavy  duty  vehicle  driving  cycle 
and  for  BMW  passenger  car  driving  cycle.  The  simulation  results 
show  that  the  hotspot  is  located  in  the  middle  of  superca¬ 
pacitors  module.  The  cooling  system  of  supercapacitors  module 
depends  on  the  power  profile  during  the  driving  cycles,  on  the 
ambient  temperature  and  on  the  convective  heat  transfer  coeffi¬ 
cient. 
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